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O ,■ Abstract 

D | Recent results from CDF Collaboration favor a large CP asymmetry in B s — > K ix 



Jl 



decay, while the Standard Model prediction is very small. Moreover, the measurement 
of its branching ratio is lower than the Standard Model prediction based on the QCD 
factorization. We compute the gluino-mediated supersymmetry contributions to B s — > 
K^*' tt + ,K(*> p + decays in the frame of the mass insertion method, and find that for 
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C"-— ■ — I < 2, the theoretical predictions including the LR and RL mass insertion contribu- 

tions are compatible with the measurements of B s — > K~ir + decay and B° — B° mixing 
within 2a ranges. Using the constrained LR and RL mass insertion parameter spaces, 

k> , we explore the supersymmetry mass insertion effects on the branching ratios, the direct 

CP asymmetries and the polarization fractions in B s — > K*~tt + , K~p + ,K*~p + decays. 
We find the constrained LR and RL insertions can provide sizable contributions to the 
branching ratios of B s — > K*~tt + ,K^*' p + as well as the direct CP asymmetry and the 
longitudinal polarization of B s — > K*~ p + decay without conflict with all related data 
within 2cr ranges. Near future experiments at Fermi Lab and CERN LHC-b can test our 
predictions and shrink/reveal the mass insertion parameter spaces. 
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1 Introduction 



In the recent ten years, the successful running of B factories BABAR and Belle has provided rich 
experimental data for B ± and B°, which has confirmed the Kobayashi-Maskawa CP asymmetry 
mechanism in the Standard Model (SM) and also shown hints for new physics (NP). Among the 
rich phenomena of B decays, the decay modes of B mesons into pairs of charmless mesons are the 
known effective probes of the CP violation in the SM and are sensitive to potential NP scenarios 
beyond the SM. The two body charmless B s decays will play the similar role in studying the 
CP asymmetries (CPA), determining CKM matrix elements and constraining/searching for the 
indirect effects of various NP scenarios. Recently the CDF Collaboration at Fermilab Tevatron 
has made the first measurement of charmless two-body B s — > K~-k + decay |lH4| 



B(B S -»■ K~7r + ) = (5.0 ± 0.7 ± 0.8) x 10" 
A d j P (B s -4 K~ir + ) = 0.39 ± 0.15 ± 0.08. 



(1) 



The measurement is important for understanding B s physics, and also implies that many B s 
decay modes could be precisely measured at the LHC-b. 



Compared with the theoretical predictions for these quantities in Refs. 



the QCD factorization (QCDF) [8|, the perturbative QCD (PQCD) [9J, and the soft-collinear 
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5j-[7|, based on 



effective theory (SCET) 10], respectively, one would find the experimental measurement of this 
branching ratio agrees with the SM predictions with SCET [5|, but lower than the predictions 
with QCDF and PQCD |d, ItJ. For the CDF measurement of A d j P (B s ->• K~7i + ), its central 
value favors a large CP violation in B s — > K~7T + decay (different from at 2.3cr), although it 



is also compatible with zero. In Refs. 



11 



12], a robust test of the SM or a probe of NP is 



suggested by comparison of the direct CP asymmetry in B s —±K 7r + decay. 

The decays B s — > K^~n + , K^~ p + have been extensively studied in the literatures (for 
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-|l7]). The tree-dominated decays B s — > K^*' tt + , K^ p + are induced 



example, Refs. pTZI, 

by b — y uud transition at the quark level, where the direct CPA are expected to be small in 
the SM. At present, among many measurements in the similar modes of B u d decays, several 
discrepancies with the SM predictions have appeared in tree-dominated b — > uuq (q = s, d) 



processes, for example, B — > 7i7r,iiK puzzles 18I-I22]. Although the discrepancies are not 



statistically significant, there is an unifying similarity pointing to NP (for example, Refs. 23 



25] ). There could be also potential NP contributions in B s — > K^*'~tt + , K^~ p + decays. The 
measurement given in Eq. flTJ will afford an opportunity to search/constrain NP scenarios 
beyond the SM. 

Supersymmetry (SUSY) is an extension of the SM which emerges as one of the most promis- 
ing candidates for NP beyond the SM. In SUSY, supersymmetric version of the SM contributes 
to the Flavor Change Natural Current (FCNC) processes. The flavor-changing in these pro- 
cesses is intrinsically tied to usual CKM-induced flavor-changing of the SM (If that were the 
only new source of flavor physics, we would say the model is minimally flavor violating). But 
general SUSY is not minimally flavor violation. For general SUSY, a new source of flavor vi- 
olation is introduced by the squark mass matrices, which usually can not be diagonalized on 
the same basis as the quark mass matrices. This means gluinos (and other gaugios) will have 
flavor- changing couplings to quarks and squarks, which implies the FCNCs are mediated by 
gluinos and thus have strong interaction strength. In order to analyze the phenomenology of 
non-minimally flavor violating interactions in general SUSY framework, it is helpful to rotate 
the effects so that they occur in squark propagators rather than in couplings, and to parameter- 
ize them in terms of dimensionless mass insertion (MI) parameters (S^ B )ij with (A, B) = (L, R) 
and (i,j = 1,2,3). In this paper, we work in the usual MI approximation 26|, |27|, and con- 
sider B s — V K^~tt + , K(*'~ p + decays, in general SUSY models, where flavor violation due to 
the gluino mediation can be important. The chargino-stop and the charged Higgs-top loop 
contributions are parametrically suppressed relative to the gluino contributions, and thus are 



ignored following [27H30J . In our work, we also discuss the implications of B° — B° mixing since 
the relevant MI parameters, that affect B s — > K^*^tt + , K^~ p + decays, enter also in B° — B° 

o 

7Tl~ 

mixing. We consider the LR, RL, LL and RR four kinds of the Mis. We find that for — § < 2, 

1 

our predictions including the LR or RL MI effects are compatible with the measurements of 
B s — > K~tt + decay and B° — B° mixing within 2<r ranges, and the constrained both LR and 
RL Mis could significantly affect the polarization fractions of B s — > K*~ p + decay. While the 
constrained LL and RR insertions from B° — B° mixing can not explain the possible large CP 
asymmetry and the small branching ratio of B s —$■ K~ii + because of lacking the gluino mass en- 
hancement in the decay. Therefore, with the ongoing 5-physics at Tevatron, in particular with 
the onset of the LHC-b experiment, we expect a wealth of B s decay data and measurements of 



these observables could restrict or reveal the parameter spaces of the LR and RL insertions in 
the near future. 

The paper is arranged as follows. In Sec. 2, the relevant formulas for B s — > K^*^ii + , 
K^~ p + decays and B° — B° mixing are presented. We also tabulate the theoretical inputs 
in this section. Sec. 3 deals with the numerical results. Using our constrained MI parameter 
spaces from B s — > K~n + decay and B° — B° mixing, we explore the MI effects on the other 
observable quantities, which have not been measured yet in B s — > K^~7i + , K^~ p + decays. 
Sec. 4 contains our summary and conclusion. 

2 The theoretical frame 

2.1 The decay amplitudes for B s — V K^~7r + , K^~p + decays 
2.1.1 The decay amplitudes in the SM 

In the SM, the low energy effective Hamiltonian for the b — > uud transition at the scale p ~ nib 



is given by [31] 

H s e %(AB = 1) = % £ X p (cf M Ql + C s 2 M Ql + £ C? M Q, + C™Q 7l + C^Q Sg ) +h.c, (2) 

V ^ p=u,c \ i=3 / 

here X p = V p bV* d with p g {u, c} are CKM factors, the Wilson coefficients within the SM Cf M 
can be found in Ref. 3jJ, and the relevant operators Qi are given as 



Q\ = {Pa^Lb^idpJuLpp), Ql = (Pal' I Lbp)(d l3 'y l ,Lp a ), 

Qs = (da'fLba) Y,%lMp)i Q* = {dprfLba) Y^cHu 1 ^^ 

q* q* 

q> q> 

3 3 - 

Qi = 7,(d a ^Lb a )Y,e q '(qpl»R(l'(3)i Qs = idd^Lbaj^eq.^^Rq'p), 

q' q' 

3 3 - 

Qg = ^(da^Lb^Y.eq'iq^pLq^), Q w = -(dpfLba) J2 e AfcyMp)> 

q' q' 

Qr-y = ^rn b d a a^Rb a F^, Q 8g = ^m b d a a^RT^bpG%, (3) 

where a and (3 are color indices, and L(R) = (1 — (+)75). 



With the weak effective Hamiltonian given in Eq. (J2J), one can write the decay amplitudes 
for the relevant two-body hadronic B —} M\M 2 decays as 

A SM (B ^M l M 2 ) = ( x M l M 2 \H^j(AB = 1)\B] 

= Y,EKCr i ^)(M 1 M 2 \Q^)\B). (4) 

P i 

The essential theoretical difficulty for obtaining the decay amplitude arises from the evaluation 
of hadronic matrix elements {MiM 2 \Qi{jj)\B) , for which we will employ the QCDF [§| through- 
out this paper. We will use the QCDF amplitudes of these decays derived in the comprehensive 



papers [6, 



13| as inputs for the SM amplitudes. 



2.1.2 SUSY effects in the decays 



In SUSY extension of the SM with conserved R-parity, the potentially most important contribu- 
tions to Wilson coefficients of penguins in the effective Hamiltonian arise from strong- interaction 
penguin and box diagrams with gluino-squark loops. They can contribute to the FCNC pro- 
cesses because the gluinos have flavor- changing coupling to the quark and squark eigenstates. 
In general SUSY, we only consider these potentially large gluino box and penguin contribu- 



tions and neglect a multitude o 



other diagrams, which are parametrically suppressed by small 



electroweak gauge coupling 27H30|. The relevant Wilson coefficients of b — ¥ uud process due to 



the gluino box or p eng uin diagram involving the LL and LR insertions are given (at the scale 



ji ~ mw ~ iTiq) by [27 
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SUSY, 



rrin 



rrin 



rriq) 



a 2 s {m q ) 



2y/2G F \ t m 2 d V 9 



■-B 1 (x) - -B 2 (x) - —Pi(x) 
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2^00 ) (5 a LL ) 13 , 



a 2 s (rriq) 



2V2G F X t m 2 d 



-\bi{?) 



l -B 2 {x) + ip^s) + \p 2 {x) ) (<^) 13 , 



a 2 s (rriq) 
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2y/2G F X t m 2 d V 9 



B x {x) + -B 2 (x) - -P^x) - -P 2 (x)) (S d L ) l3: 



a 



SUSY I 



a 2 s (mq) 



2 
2V2G F \ t mj, \~3 Bl ^ X ' 



6 



~*SUSY ( 
'l-i \" b q) 



Cv~~- (rrin) 



a 



SUSY, 



m h 



9V2G F X t m1 
2-Ka s (mg) 
V2G F X t m 2 d 



B 2 (x) + -PxOc) + -P 2 {x) J (S d LL 



(13) 



(^) 13 M 4 (x)-(^) 13 ^j4S 1 (x 



(6 d LL ) 13 (-M 3 (x) - -M A (x 

H5 d LR ) 13 Q-^(AB 1 (x)-9x~ 1 B 2 (x: 



(5) 



where x = rn?~/m 2 ~, and the loop functions B^x), Pt(x), M^x) can be found in Ref. 32[. For the 



RR and RL insertions, we have additional operators Qi=3... 6,77,83 that are obtained by L -h- i? 
in the SM operators given in Eq. (J2J). The associated Wilson coefficients C^ 3 sl ^ 7l8g are 
either dominated by their expressions as above with the replacement L -h- R. The remaining 
coefficients are either dominated by their SM (C\^) or are electroweak penguins (C7...10) and 
therefore small. 



The SUSY Wilson coefficients at low energy Cf (// ~ m&) can be obtained from C, 



iSUSY, 



rrin 



in Eq. (J5J) by using the renormalization group equation as discussed in Ref. [31 1 



C(fi) = U 5 (fi,mg)C(mq) 



qji 



(6) 



where C is the 6x1 column vector of the Wilson coefficients and [^(/x, m g -) is the five-flavor 6x6 
evolution matrix. The detailed explicitness of U^(u, mq) is given in Ref. 3jJ. The coefficients 
Cj} JSY and C , fj 7sy at the /j ~ ra& scale are given by 
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riSUSY 

°7 7 

riSusY 

°8g 



(A*) 
(A*) 



where 77 



'°>N^f».M^ 



■a«W ' 



*a„( m b)* 






(7) 



2.1.3 The total decay amplitudes 

For the LL and LR insertion, the NP effective operators have the same chirality wit 
the SM, so the total decays amplitudes can be obtained from the SM ones in Refs. 
replacing 



1 ones of 



« 



Uby 



ct M ->• ct M + a 



SM 



iSUSY 



(8) 



For RL and RR insertion, the NP effective operators have the opposite chirality with the SM 
ones, and we can get the corresponding decay amplitudes from the SM decay amplitudes by 



following replacements 



37] 



ct M -> a 



SM 



a 



SUSY 



(9) 



for A(B S ->• K~tt + ) and A 0: \\(B S -> K*~p + ), as well as 



iSUSY 



cr -+ Cf M + Ct 



(10) 



for A(B S -»■ K*-tt + ), A(B S -> A'~p + ) and A ± (5 S -»■ AT*V 
Then the total branching ratios read 

Tb, \p. 



B(B S -> MiM 2 ) = -^- |^(S a -> M!M 2 )| Z , (11) 



57Tm 



where tb s is the B s lifetime, \p c \ is the center of mass momentum in the center of mass frame 
of B s meson. 

The direct CP asymmetry is defined as 

di B(B S -> /) - B(B S -> /) 

^ p(jBs ~* /} " B(B a ^f) + B(B 8 ^fY (12) 

In the £> s — > yy decay, the two vector mesons have the same helicity, therefore three 
different polarization states are possible, one longitudinal and two transverse, and we define 
the corresponding amplitudes as Ao t ±. Transverse («4.||,j_) and helicity (A±) amplitudes are 
related by A\\ t ± = + y=? - . Then we have 

\a{b s ^w)\ 2 = iAi 2 +iAr+i^-r = iAi 2 +i^ni 2 +i^±r. (13) 

The longitudinal(transverse) polarization fraction fi(f±) are defined by 



IV IAj i! 

\Aq\ 2 + \A l{ \ 2 + \A A 



hABs-^w) -^= lAl7 .'r:,, 1 , ,, , 2 - (14) 



2.2 Effective Hamiltonian for B° — B° mixing 

The most general B° — B° mixing is described by the effective Hamiltonian 

H eff (AB = 2) = J2 C[Q\ + £ C[Q' t + h.c, (15) 



j=l i=l 



with 



Qi = (d-fPLbMdjpPLb)!, 

Q' 2 = {dP L b)x{dP L b)i, 

Q' 3 = (dP L b) s (dP L b) 8l 

Q* = {dPLb)i{dP R b) x , 

Q' 5 = (dP L b) 8 (dP R b) 8 , (16) 



where Pl{r) = (1 _ (+)7s)/2 and the operators Q\ 2 3 are obtained from Q[ 2 3 by the exchange 
L <r+ R. The hadronic matrix elements, taking into account for renormalization effects, are 
defined as 



with S Bd = ( 



(B°\Q'M\B°) = 


^n 2 B jl d B x {n), 


(B°\Q'M\B°) = 


5 

-J^™l d fB d S B d Mri 


(B°\Q'M\B°) = 


Y^ m B d fB d SB d B 3 (fl), 


(B°\QM\B°) = 


2 m B d fB d SB d B4fj,), 


(B°\Q' 5 (ri\B°) = 


Q m B d fB d S B d B 5 (fJ,), 


m B d \ 2 





(17) 



3 d \m b (m b )+m d (m b ) J 

The Wilson coefficients C[ receive contributions from both the SM and the SUSY loops: 
C[ = C'f^ 1 + C'f . In the SM, the t — W box diagram generates only contribution to the 
operator Q[, and the corresponding Wilson coefficient C[ SM at the m b scale is 31 ] 



G 



C[ SM (m b ) = -^m 2 w (V td V t l) 2 V2B S (x t )[a s (m b )r^ 



471-2 



a s (m b ) 



(18) 



where x t = mf/m^ and r]2B is the QCD correction. 

The gluino-mediated SUSY contributions to B°—B® mixing and 2/3 in the MI approximation 
have been extensively studied (for example, Refs. 



38 



40- 



42]). In general SUSY models, 
there are new contributions to B° — B° mixing from the gluino-squark b ox diagrams, and the 



corresponding Wilson coefficients C' i SUSY (at the rriq scale) are given by 27H30] 



cf^K~) 



C' 2 SUSY (m,) 

,r K) 



o-r 



216m? 



o; 



216m? 



(24x/ 6 (x) + 66/ 6 (x))(^ L )? 3 , 



204xf 6 (x)(8 d RL )l 3 , 



-if SUSY i 
'3 



o; 



216m? 



36xf 6 (x)(S d RL ) 2 X3 , 



Ci SUSY (m,) 



C' 5 SUSY (m,) 



ri 



216m? 



(504x/ 6 (x) - 72/ 6 (x)) (^)i 3 (4h)i3 - l32f 6 (x)(5 d LR )i3(S d RL ] 



13 



n: 



216m? 



24x/ 6 (x) + 120/ 6 (x) (^ l ) 13 (4h)i3 - 180; 6 (*)(<5£ r )i3(4l)i3 (19) 



The loop functions f&(x), fe(x) can be found in Ref. 32]. The other Wilson coefficients C[ S 2 U 3 SY 
are obtained from C[ 23 by exchange oi L ^ R. 



The SUSY Wilson coefficients at the mj, scale Cf USY (mb) can be obtained by 

C r (m b ) = E E (t s) + V'cp s) ) fC s K), (20) 



where r/' = a s {mq)/a s {m t ). The magic number aj b[ and c[ can be found in Ref. 381 ]. 
Renormalization group evolution of Ci.2,3 is done in the same way as for Ci,2,3- 

In terms of the effective Hamiltonian in Eq. (J15I) . the mixing amplitude M12 read 



= W«.//(A^2)|^) 

2 ™B d 

Then, the _B° mass difference AM^ = 2|Mi2|, its associated CP phase 2(3 = arg(Mi2), and 
S^k s = sin2/3. Experimental values of AM^ and 2/3 are given by the Heavy Flavor Averaging 



Group [39 ] 



AM d = 0.507 ± 0.005 ps _1 , 

sin2/3 = 0.68 ±0.03. (22) 

Above experimental bounds will also be used in our numerical results. 

2.3 Input Parameters 

The input parameters are collected in Table [TJ We have several remarks on the input parame- 
ters: 

• Wilson coefficients : The SM Wilson coefficients Cf are obtained from the expressions 
in Ref. [31 1 . 

• CKM matrix element : For the SM predictions, we use the CKM matrix elements from 



the Wolfenstein parameters of the latest analysis within the SM in Ref. 44j, and for 



the SUSY predictions, we take the CKM matrix elements in terms of the Wolfenstein 



parameters of the NP generalized analysis results in Ref. 44]. 



Masses of SUSY particles: When we study the SUSY effects, we will consider each possi- 
ble MI (<^4#)i3 for AB = LL, LR, RL, RR only one at a time, neglecting the interferences 
between different insertions products, but keeping their interferences with the SM am- 
plitude. We fix the common squark masses nig = 500 GeV and consider three values of 
x = 0.25, 1,4 (i.e. m- g = 250, 500, 1000 GeV) in all case. 



Table 1: Default values of the input parameters. 



m w = 80.398 ±0.025 GeV, m Bd = 5.280 GeV, m Bs = 5.366 GeV, r Bs = (1.472±g;gl|) ps, 

m K , ± = 0.892 GeV, m i<± = 0.494 GeV, m^ ± = O.WOGeV, m p = 0.775 GeV, 

m t = 171.3±?;| GeV, m b (m h ) = (4.20 ± 0.07) GeV, m„(2GcV) = (0.0015-0.0033) GeV, 

m d (2GeV) = (0.0035 - 0.0060) GeV, m s (2GeV) = (0.105+g;^) GeV. [43J 

The Wolfcnstein parameters for the SM predictions: 

A = 0.810 ± 0.013, A = 0.2259 ± 0.0016, p = 0.154 ± 0.022, 77 = 0.342 ± 0.014. 

The Wolfenstein parameters for the SUSY predictions: 

A = 0.810 ± 0.013, A = 0.2259 ± 0.0016, p = 0.177 ± 0.044, fj = 0.360 ± 0.031. [44] 

f K = 0.160 GeV, f K , = (0.217 ±0.005) GeV, f£ t = (0.156 ± 0.010) GeV, 

f„ = 0.131 GeV, f p = (0.205 ± 0.009) GeV, ff = (0.147 ± 0.010) GeV, 

A B S ^K* ^ = q 36Q ± q Q34) A B,^K* ^ = q 233 ± q Q 22) A B S ^K' ^ = Q lgl ± q Q25) 

V B '^ K '(0) = 0.311 ±0.026, f£--* k (Q) = 0.301^3. [45, 46] 



f Bs = (0.245 ±0.025) GeV, f Bd = (0.200 ± 0.020) GeV, f Bd \J B Bd =0.225 ±0.025 GeV. 


[47] 


X B = (0.46 ±0.11) GeV. 


[48] 


r] 2B =0.55 ±0.01. 


[49] 


a \ =0, al = 0.20 ± 0.15, al = 0, a p 2 = 0.1 ± 0.2, 

af = 0.2 ± 0.2, of = 0.1 ± 0.3, af * = 0.06 ± 0.06, af * = 0.1 ± 0.2. 


[6,i3] 


B<f\m b )= 0.87(4) (j*) , B^(m 6 ) = 0.82(3)(4), B^(m 6 ) = 1.02(6)(9) 
Bf (m b ) = 1.16(3) (t 5 7 ) , Bf(m 6 ) = 1.91(4) (+») 


[50] 



3 Numerical results and analysis 

Now we are ready to present our numerical results and analysis. First, we will show our 
estimations in the SM with the parameters listed in Table [T] and compare with the relevant 
experimental data. Then we will study the SUSY predictions for the branching ratios, the CP 
asymmetries and the polarization fractions in B s — > K^~n + , K^~ p + decays. For CPA of 
B s — > K*~ p + , we will only study the longitudinal direct CP asymmetry (A c ' P ir ). 

In the SM, the numerical results with la error ranges for the sensitive parameters are 
presented in Table [2j The detailed error estimates corresponding to the different types of 
theoretical uncertainties have been already studied in Refs. la [131, and our SM results of B, 

fin 

A d Q P and fi are consistent with the ones in Refs. [y, |13|]. For the color-allowed tree-dominated 
decays B s — > K~ n + ,K*~ n + ,K~ p + ,K*~ p + , power corrections have limited impact, and the 

10 



Table 2: The SM predictions with ler error ranges of the input parameters for B (in units of 
KT 6 ), Af P (in units of KT 2 ) and f L in B s ->■ AW-vr+, K^~p+ decays within QCDF. 



Decay modes 


B 


Adir ( AL,dir\ 

•^■cpy ^CP J 


h 


B s -> A^tt+ 


[6.89, 15.67] 


[-8.11,-1.44] 




5 S -»■ /i*-7T + 


[9.03, 18.36] 


[-0.64,2.01] 




£? s -)• Ji- P + 


[14.79,38.39] 


[-1.68,0.86] 




£? s -»• A— p+ 


[11.85,69.69] 


[-5.19,-1.02] 


[0.87, 0.97] 



main sources of theoretical uncertainties in the branching ratios are the CKM matrix elements 
and the form factors. Their Aq P and A(jp r can be predicted quite precisely, and found to be 
very small (~ 10~ 2 ) due to small penguin amplitudes. The uncertainty of Jl{B s — > K*~ p + ) 
is mostly due to the uncertainties of different form factors. Comparing the SM predictions in 
Table [2] with the relevant experimental data in Eq. (TO , we may see the present experimental 
data of the B s — > K~ix + decay within ler ranges are not consistent with the SM predictions 
within ler error ranges of the input parameters. The experimental data of Aq P (B s —$■ K~n + ) 
are obviously larger than the SM prediction, moreover, the measurement of its branching ratio 
is lower than the SM prediction based on the QCDF. 

Now we turn to the gluino-mediated SUSY contributions to B s — > K^~tt + , K^~ p + decays 
in the framework of the MI approximation. These decays are induced by b — > uud transition 
at the quark level, and we consider four kinds of Mis (LL, LR, RL and RR) contributing to 
B s — y A(*)~7r + , K^*'~ p + decays. In the SM, the very small direct CPA of these decays come 
from the weak phase of small penguin amplitudes. In order to have nonzero CPA, we need at 
least two independent amplitudes with different weak phases. In general SUSY models, we are 
considering, the weak phases reside in the complex MI parameters Ss and appear in the SUSY 
Wilson coefficients in Eq. (jSJ). These weak phases are odd under a CP transformation. 

First, we study constraints on the MI parameters (Sf L )i3, (5f R )i3, (S RR )i3 and (6 RL )i3 
from B s —y K~tt + decay and B° — B° mixing. Using the formulas in Sec. 12.11 and the input 
parameters in Sec. 12.31 we may obtain the two-body branching ratios and the direct CPA 
of B s — >• K ( -*'~7i + , K^*'~ p + decays within the framework of SUSY by the QCDF. The mass 



11 



difference and CP phase of B° — B° mixing are gotten from the formulas in Sec. 12.21 and the 



input parameters in Sec. 12.31 Noted that we take the relevant CKM matrix elements in terms of 



the Wolfenstein parameters of the NP generalized analysis results in [44j, which is different from 
the SM predictions. In each of the MI scenarios to be discussed, we vary the mass insertions 
over the range |(^xs)i3| < 1 to fully map the parameter space. We then impose experimental 
constraints from B s — > K~ir + decay and B° — B° mixing, which are shown in Eq. (FjQ) and Eq. 
respectively. 
For the LL and RR Mis, |(<5^ L )i3| and | (^^j?)i3 1 are strongly constrained from B° — B° 
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38 



40 



41] . The effects of the constrained LL and RR insertions on B s — > K~ix + 
are almost negligible because of lacking the gluino mass enhancement in the decay, and they 
will not provide any significant effect on B(B S — > K~tt + ) and Aq P (B s -» K~ir + ). The current 
data within lcr (or 2<r) ranges of Aq P (B s — >• K~n + ) cannot be explained by the LL and RR 
Mis. 

The case of the LR or RL insertion is very different from that of either LL or RR insertion. 
The LR and RL Mis only generate (chromo)magnetic operators Qj^fig an d Q7 7 ,8g 5 respectively. 
Especially, the LR and RL insertions are more strongly constrained, since their contributions 
are enhanced by nig/mb due to the chirality flip from the gluino in the loop compared with the 
contribution including the SM one. In these cases, even a small (c^r)i3 or (S PL )i 3 can have 
large effects in the decays. 

We can not get the allowed spaces of the LR and RL insertion parameters which may explain 
all these data within la ranges given in Eq. (TjQ) and Eq. f[2"2"j) at the same time. In fact, there 
is no common allowed parameter space from B(B S — > K~n + ) and Aq t p (B s — > K~n + ) within 
lcr ranges. Using the central values of the input parameters and nig = m g = 500 GeV, we 
show the contour plots of B(B S — > K~n + ) and Ac P (B s — > K~tt + ) in the K^^)^! — 4>lr and 
I(^l)i3| — ^RL planes in Fig. [TJ We discuss Fig. d](a) in detail. Fig. d](a) show the contour plot 
of la error bars of B(B S — >• K~n + ) and Aq P (B s — > K~n + ) in | (5^)13 1 — 4>lr plane. The space 
between the contour lines of B(B S ->■ K~ix + ) = 6.1 x 10~ 6 and B(B S -> K~n + ) = 3.9 x 10" 6 
is the allowed space for B(B S — > K~n + ) within lcr ranges. The space between the contour 
lines of AfJ P (B s -> A'-vr + ) = 0.56 and A d j P (B s -> A'-vr + ) = 0.22 is the allowed space for 
Aq P (B s — y K~tt + ) within lcr ranges. From Fig. Q](a), we can see that there is no intersection 
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Figure 1: The contour plots of i3(5 s ->■ #~7r+) and A% r P {B s -> #-7r + ) in the |(c^) 13 | - Lfl 



and |(4 L ; 



13 I 



!»_rl planes with the central values of the input parameters and rriq 



m- n 



500 



GeV. Plots (a) and (c) show the la error-bar of data given in Eq. (JTJ). Plots (b) and (d) show 
the 2a error-bar of data given in Eq. (TO . 0^ denotes the mixing parameters weak phase. 
Please see text for details. 

between the allowed space from B(B S — > K~tt + ) and the allowed space from Aq P (B s — > K~tt + ), 
and the experimental data of A^ r P (B s — > K~n + ) within lcr give quite strong constraints on 
(^Li?,)i3- The similar results for the RL insertion are shown in Fig. [1] (c) in 1(5^) 13] — 4>rl 
plane, and there is also no intersection between the allowed space from B(B S — > K~ir + ) within 
lcr ranges and the allowed space from Aq P (B s — > K~ti + ) within la ranges. 

Then we expand the experimental bounds within 2a ranges to search for the allowed spaces 
of the LR and RL MI parameters. As shown in the dark gray ranges of Fig. [1] (b) and (d), there 
are the allowed spaces for the LR and RL insertion parameters from the data of B s — > K~tt + 
decay within 2a ranges. Noted that the allowed parameter spaces shown in the dark gray 
ranges of Fig. [1] (b) and (d) are obtained by using the central values of the input parameters. 
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The allowed spaces will be enlarged if we consider the theoretical uncertainties of the input 
parameters. We use the experimental bounds of B s — Y K~ir + decay and B° — B° mixing within 
2a ranges, and take the input parameters within la ranges to obtain the allowed spaces of the 
LR and RL insertion parameters. The constrained spaces of (cr| i? ) 13 and (cr^ L ) 13 for rrig = 500 
GeV and different x are demonstrated in Fig. [21 and the corresponding numerical ranges are 
summarized in Table [31 
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Figure 2: The allowed parameter spaces of the LR and RL MI parameters constrained from 
B s —?• K~ir + decay and B° — B° mixing at 95% C.L. for the squark mass rrig = 500 GeV and 
for the different values of x = 0.25 and x — 1. 



Table 3: Bounds on the LR and RL MI parameters from the measurements of B s — > K tt + 
decay and B° — B° mixing at 95% C.L.. 



x 



0.25 
1 



I(«£h)is| 



4>lr (deg.) 



[0.005,0.031] 
[0.008, 0.056] 



-85,15] 
-81,10] 



m 



RLm\ 



<f>RL(deg. 



[0.005, 0.030] [100, 180] U [-180, -158] 
[0.008, 0.054] [94, 180] U [- 180, - 154] 



In Fig. 121 we can see the allowed moduli of the LR and RL MI parameters are very sensitive 
to the values of x, nevertheless the allowed phase ranges of the LR and RL MI parameters are 
not changed obviously for different x. We find that both moduli and phases of {Slrrl)i3 are 
strongly constrained by the data of the B s — > K~ir + decay within 2a ranges, since the SM 
prediction of A.q P (B s — > K~n + ) with the QCDF is not consistent with the experimental data 
within 2a ranges. The lower limits of |(^l,r.rl)i3| come from the experimental lower limit of 
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Afj P (B s —?• K~ir + ) within 2a ranges. In the case of x — 1, 4>lr,rl are also constrained by 
AMd within 2a ranges, while sin2/3 within 2a ranges doesn't provide any further constraint. 
It's different in the case of x = 0.25, <Plr,rl are also constrained by sin2/3 within 2a ranges, 
but AM,j within 2a ranges doesn't provide any further constraint. 

It is worth to note that we also study the constrained spaces of (5^) 13 and (S RL )i3 in 
the case of x — 4, and we find there is no intersection of the constrained spaces between 
from B° — B° mixing within 2a ranges and from B s — > K~n + decay within 2a ranges. From 
B° - B° mixing within 2a ranges, we get \(5i RtRL )i 3 \ < 0.07, in which <p LRRL G [-180°, 180°] 
if \Kr,rl)i3\ < 0-03, and <\> LRRL G [-110°, -45°] U [70°, 135°] if \{5 d LR>RL )i 3 \ G [0.03,0.07]. We 
obtain |(<51rr£)i3| G [0.03,0.19] from B s — > K~n + decay within 2a ranges, in which <pL R G 
[-90°, 5°] ((f> RL G [-180°,-175°]U[90°,180°]) if |(^ L ) 13 | G [0.07, 0.19], and <j> LR G [-40°, 15°] 
{<f>RL G [-180°, -165°]U[140°, 180°]) if \(S d LRRL ) 13 \ G [0.03,0.07]. So there is no common allowed 
phases between the bounds from B° — B° mixing and the bounds from B s — > K~ir + decay in 
the case of x — 4. If x — 2, there still are common allowed spaces for the LR and RL Mis. 
For the LR MI, two very narrow allowed spaces are near (|(^ R )i3| ~ 0.02, <^# w 4°) and 
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(\(^lr)is\ ~ 0.04, 4>lr ~ —64°). For the RL MI, two very allowed spaces are near (\(S RL/ 
0.02,0^ « -176°) and (|(^)i 3 | w 0.04,0^ w 116°). 

The relevant upper bounds have been obtained in Refs. 3§J, |40|, |4l|. In Ref. [38j, 
\(3lrrl)i3\ — 0.07, 0.08, 0.11 for x = 0.25, 1, 4, respectively, which are constrained by 
imposing the experimental bounds from AMd and S^k s within la ranges, setting rriq = 500 



GeV and scanning over the CKM phase 7 G [0, 27r]. In Ref. |4QL,J(5£ flBL )i 3 | < 0.015 for 



m - = m- = 500 GeV from AM d and 2/3 within 2a ranges. In Ref. [4l|, \(5 RLRL )i 3 \ < 0.03 for 



nig, rriq < 600 GeV from AMd and 2/3. In the case of x = 0.25 and x = 1, comparing with the 



exist bounds in 



38 



40 



4l|, our upper limits of 1(5^^) 13 1 are at the same order of previous 
ones, while the lower limits of |(<5iij.zyji3| are also given from B s — > K~tt + decay within 2a 
ranges. In addition, 4>lr,rl are strongly constrained from B s — > K~ir + decay. The allowed 
space for x = 4 case are ruled out by both B s — > K~ir + decay and B° — B° mixing together. 

Next, we will explore the SUSY effects on the other quantities, which have not been mea- 
sured yet in B s — > K*~tt + ,K^~ p + decays, by using the constrained parameter spaces of the LR 
and LR insertions as shown in Fig. [2j With the expressions for B, Afjp and /l, we perform a 
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scan through the input parameters within 2er ranges and the new constrained SUSY MI param- 
eter spaces, and then the allowed ranges for B, Aqp and J'l are obtained with different SUSY 
mixing insertion parameter, which satisfy relevant experimental constraints of B s — > K~tt + 
decay given in Eq. ([T]) and B° — B° mixing given in Eq. (1221) . The numerical results for 
B s —?• K*~7c + , K^*'p + with different x value are summarized in Table HI The corresponding 
SM predictions with 2<r error ranges of the input parameters are also listed for comparison in 
the second column of the Table HI We can see the data of B{B S — > K~n + ) is consistent with the 
SM prediction of B(B S —$■ K~n + ) at 95% C.L., nonetheless very close the lower limit of its SM 
prediction. The data of Ac P (B s — > K~ti + ) is not consistent with its SM prediction at 95% C.L.. 
From the last four columns of the Table HI we can see the results are similar for different value 
of x. In the SUSY predictions of the branching ratios, there are many sources of uncertainties, 
mainly arising from different form factors, CKM matrix elements, the annihilation contribution, 
other hadronic parameters of the QCDF, and the constrained MI parameters. The uncertain- 
Table 4: The theoretical predictions for B (in units of 10~ 6 ), Aq P (in units of 10 -2 ) and /x, 
in four B s — > K^~ir + , K^*'~ p + decays based on general SUSY models with different MI and 
different x. The corresponding SM predictions with la error ranges of the input parameters 
are also listed for comparison. 



Observables 


SM predictions 


SUSY values 

with (Sf R )i3 

for x = 1 


SUSY values 
with (6i R )i 3 
for x = 0.25 


SUSY values 

with (<5^, L )i3 

for x = 1 


SUSY values 
with (S RL )i3 
for x = 0.25 


B(B S -> K-7T+) 


[3.32, 24.48] 


[2.80,7.20] 


[2.80,7.20] 


[2.80,7.20] 


[2.80,7.20] 


B(B S -> A"*-7T+) 


[5.91,25.80] 


[0.25,30.67] 


[0.13,30.94] 


[7.10,93.54] 


[7.53,93.32] 


B(B S -> K~p+) 


[7.36, 57.97] 


[0.01,33.31] 


[0.01,31.51] 


[15.96,275.11] 


[17.51,273.01] 


B{B S -> K*~p+) 


[2.47,127.56] 


[1.36,64.77] 


[1.21,77.47] 


[1.21,72.76] 


[1.08,81.71] 


A C p(B s -^K-tt+) 


[-14.46,0.97] 


[5.00,25.07] 


[5.00,22.02] 


[5.00,22.73] 


[5.00,19.10] 


A C p(B s -»■ K*-n+) 


[-4.62,6.12] 


[-12.31,4.88] 


[-20.27,13.30] 


[-5.15,4.84] 


[-4.40,5.21] 


A C p(B s -> K-p+) 


[-3.63,2.03] 


[-36.59,13.70] 


[-35.67,34.38] 


[-1.85,1.48] 


[-1.76,1.26] 


A L c f\B s -> K*-p+) 


[-13.32,1.62] 


[-80.64,36.23] 


[-94.55,63.70] 


[-54.42,56.75] 


[-75.29,65.98] 


f L (B s ^K*-p+) 


[0.51,0.98] 


[0.02,0.97] 


[0.01,0.98] 


[0.03,0.98] 


[0.01,0.98] 
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ties of these direct CPA mostly come from the constrained MI parameters. The uncertainty of 
fi{B s — > K*~ p + ) due to different form factors and the constrained MI parameters. 

Comparing the SUSY predictions to the SM predictions given in Table HJ we give some 
remarks on the numerical results: 

• The LR and RL Mis have significant effects on B(B S — > K~tt + ), and the relevant pa- 
rameters have been limited by both upper and lower limits of B(B S — > K~n + ) within 2a 
ranges. The LR and RL Mis also have great effects on A^ P (B S — > K~tt + ), which could 
be increased from the SM prediction range [—0.14, 0.01] to about [0.05, 0.2], however, the 
range [0.05, 0.2] is still far from the central value of its measurement, and is near to the 
lower limit of the measurement within 2a ranges. The LR and RL MI parameters just 
have been limited by the lower limit of A^ P (B S — > K~ir + ) within 2a ranges. 

• The constrained LR insertion still has significant effects on B(B S — > K*~ir + , K^~ p + ). 
The allowed lower limits of B(B S — > K*~tt + , K~ p + ) could be reduced one or two order(s) 
from their SM predictions, and the allowed upper limit of B(B S — > K*~ p + ) has been 
suppressed a lot from its SM prediction. In the SM, the direct CPA are very small in 
these decays. It is interesting to find the contributions of the constrained LR insertion 
have great effects on Ac P (B s — > K*~n + , K ( -*^~p + ). The allowed ranges of Ac P {B s — > 
K*~n + , K^~ p + ) could be extremely enlarged from ones of their tiny SM predictions. We 
find the constrained LR insertion contributions have a great impact on the longitudinal 
polarization fraction /l(5 s — > K*~p + ), which could be reduced to about zero by the 
constrained LR insertion. 

• For the RL MI case, the effects of the constrained RL insertion could exceedingly increase 
the allowed upper limits of B(B S — > K*~7i + ,K~p + ) and decrease the upper limit of 
B(B S ->• K*~p + ). The RL insertion has small effects on Af, r P (B s ->• A'*~7r+, K~p + ), 
however, this insertion can greatly affect A CP ir (B s — >■ K*~ p + ). /l(B s —> K*~ p + ) also 
could be reduced to about zero by the constrained RL insertion. 

Noted that the LR and RL Mis only generate dipole operators Qi-yfig and Qi^^g, respectively, 
and Q7 1 ^g ) Qi 1 ^g do not contribute to the transverse p eng uin amplitudes at 0(a s ) due to 



angular momentum conservation in B s — > K* p + decay 5JJ. In other words, the LR and RL 
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Mis only contribute to the longitudinal penguin amplitude at 0(a s ). Because the LR and RL 
contributions are enhanced by m g /mb, even a small (<^r)i3 or (5^ L )i3 can have large effects on 
the longitudinal penguin amplitude, and then can significantly affect the polarization fractions 
of B s — y K*~ p + decay. For the similar reason, the LR and RL Mis have been proposed as a 
possible resolution to the polarization puzzle in B — y <pK* decays 52|, 153 ]. 

For each LR and RL insertions, we can present the distributions and correlations of B, Aqp, 
fi within the modulus or weak phase of the constrained MI parameter space in Fig. |5]by two- 
dimensional scatter plots. The LR MI effects on all observables of B s — y K^~ir + , K^~p + are 
displayed in Fig. |3]and Fig. HI Fig. |3] shows the sensitivities of all observables to |(<^r)i3| for 
different value of x, and we see that all observables expect B(B S — y K~ir + ) are a little sensitive 
to both |(^ R )i 3 | and the values of x. Fig. @] displays the sensitivities of the observables to (J>lr 
for different x, and we can see that the weak phase (J)lr for different x value has similar allowed 
ranges, and has similar effects on every observable. In addition, for comparing conveniently, 
we show the SM bounds of these observables by orange horizontal dash lines and the limits 
of the measurements of B s —y K~n + within la error-bar by the cyan horizontal solid lines. 
From Fig. [3]^a-d) and Fig. H^a-d), we see that B[B S —y K~tt + ) is strongly constrained from its 
experimental data, B(B S — y K*~ ~tt + , K^*>~ p + ) are very sensitive to both |(5£r)i3| and 4>lr, and 
they are decreasing with \(Sf R )i 3 \ but increasing with 4>lr- As shown in Fig. [3](e) and Fig. IH(e), 
the LR insertion has positive effects on A^ r P (B s —y K~n + ), and there is no any point in the SM 
area since AQ r P (B s —y K~ir + ) is strongly constrained by the corresponding experimental data, 
which are not consistent with the SM predictions at 95% C.L.. Fig. [3] (f-h) and Fig. H] (f-h) 
display that Aq p (B s —y K*~n + , K^~ p + ) are sensitive to |(^ R )i 3 |, and could have very large 
allowed ranges when 4>l P e [—80°, —40°]. As for the LR insertion effects on fi(B s — y K*~ p + ), 
we show it in Fig. [3] (i) and Fig. H] (i), and we can see fh{B s —y K*~ p + ) could be hugely 
affected by the LR MI. Jl{B s — > K*~ p + ) has some sensitivities to both |(^ fi )i3| and 4>lr, an d 
it has smaller allowed range with <Plr- So the future measurement of Jl{B s —y K*~ p + ) could 
give obvious constraint on c/)lr- 

Next, we discuss the RL MI effects on all observables in B s —y K^~n + , K^~ p + decays. Fig. 
E]and Fig. [Hlshow the observables as functions of |(<5^ L )i 3 | and <Prl, respectively. Fig. E](a,d) 
and Fig. E]^a,d) show the constrained RL MI has negative effects on B(B S — y K~tt + , K*~p + ), 
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Figure 3: The effects of (S% R )is in # s -> i^W-7r + , A" ( * )_ p + decays. .4g£ and B are in units of 10 -2 
and 10 -6 , respectively. The orange horizontal dash lines denote the limits of SM predictions, and the 
cyan horizontal solid lines represent the 2a error bar of the measurements. (The same in Figs. 3KJ). 
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Figure 4: The effects of 4>lr in B s ->■ K^~ir + ,K^~ p + decays. 
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Figure 5: The effects of (5 d RL ) 13 in B s ^ K^~ n + ,K^-p + decays. 
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Figure 6: The effects of <p RL in B s -► lfW~7r + ,li'W~p + decays. 
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which is same as the LR MI effects on them. Fig. [S£b,c) and Fig. E^b,c) show us the constrained 
RL MI has very large positive effects on B(B S — > A*~7r + , K~ p + ), which is different from the LR 
ML B(B S — > A*~7r + , K^~ p + ) have some sensitivities to \{5 PL )i^\ and <Prl as displayed in Fig. 
[^b-d) and Fig. E](b-d). As shown in Fig. E]T,g) and Fig. El(f,g), the RL MI has small effects on 
Aqp(B s — > K*~n + , A'~p + ), and the RL insertion contributions can not be distinguished from 
the SM prediction. Fig. [5] (h) and Fig. [6] (h) show the RL MI effects on Af P {B s -)• A*~p+) 
are similar to the LR insertion effects on it. Afj P (B s — > K*~ p + ) could has very large allowed 
ranges when <p RL <G [100°, 140°]. As shown in Fig. El(i), /l(-B s — > K*~ p + ) could be strongly 
suppressed by the RL insertion, too. 

In addition, for the LR MI case, we show the resulting predictions for B(B S — > A~7r + ) 
vs. B(B S -> K*-7r + 1 K^-p + ) and B(B S ->■ K~ix+) vs. Afj r P (B s -> A*-vr+, A«-p+) in Fig. 
El as well as A d £ r P {B s -> R-ti + ) vs. B(B s -»■ K*~tt + , K { ^~ p + ) and A d £ r P {B s -> K~it + ) vs. 
Acp{B s -¥ K*-7T + ,K^-p + ) in Fig. [HJ In all plots in Fig. and Fig. El the black and green 
points satisfy the constraints of B s — > K~n + decay and B° — B° mixing within la ranges. 
As displayed in Fig. [7J both B(B S -> K*~ir + , K^- p + ) and Afj r P (B s -> A*-tt + , A«-p+) are 
not very sensitive to the constrained B[B S — > K~ir + ). However, as shown in Fig. [SJ both 
B(B S -»■ A'* _ 7r + , A w ~p + ) and Ac P (B s -> A*~7r + , A w ~p + ) have some sensitivities to the 
constrained ^4^p(5 s — >■ A~7r + ). Fig. [S] (a-c) show us, ,4*p(.B s — >■ A~7r + ) could has maximum 
when B(B S — > K*~n + , K^*'p + ) is near the lower limits of the SM predictions, and A d j P (B s — y 
K~7T + ) could have smaller allowed range with B(B S — > K*~tt + , A**) - p + ). Fig. [8] (d-f) indicate 
that there are some pionts accounting for large values of \Aq P (B s — > K*~ir + , K^~ p + )\ when 
Afj P (B s -»■ A-7T+) is small. 

There are similar correlations in the case of the RL MI as ones in Fig. [7] and Fig. [8] 
expect the correlations between B(B S — > K~ir + ) and Afj P (B s — > A~p + ) as well as between 
A d j P (B s -* A-7T+) and A d j P (B s -* K*~tt + , K~p + ). A d S P {B s -> K*~iv + ,K^p + ) have no any 
sensitivities to A d j P {B s -+ A~7r + ) and B{B S ->• K~tt + ). 
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Figure 7: For the LR insertion, the correlation plots between B{B S — > K 7r + ) and 
K*~ir + ,K(*' p + ) are shown in (a-c), and the correlation plots between B(B S — > K 
A% T P (B S ->• K*~-K + ,K^~p + ) are shown in (d-f). 



B(B S -> 
"7r + ) and 




10 20 30 40 

6(B s ->kV) 
'(a) 



10 20 30 40 

6(B s ->K P + ) 
(b) 



20 40 60 80 100 

B(B s ->KV) 
(c) 




-0.4 -0.2 0.0 0.2 0.4 
<( B s-> KV ) 

(d) 



-0.6 -0.3 0.0 0.3 0.6 

<(B S ->K> + ) 
(e) 



-1.0 -0.5 0.0 0.5 1.0 

A L J"(B->KV) 
(f) 



Figure 8: For the LR insertion, the correlation plots between Aq P (B s — > K ir + ) and 
K*~tt + ,K(*' p + ) are shown in (a-c), and the correlation plots between Afjp(B s — > K~ 
A d S r p{B s -)■ K*-tt + ,K^- P + ) are shown in (d-f). 
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4 Conclusions 

Motivated by recent results from CDF, which favor a possible large CP asymmetry and a small 
branching ratio in B s — > K~n + decay, we have studied the gluino-mediated SUSY contributions 
with the Mis to four B s — > K^~ 'ix + , K^*'~ p + decays based on the QCDF approach. For the 
LL and RR Mis, we have found that the constrained LL and RR insertion effects from B° — B° 
mixing are almost negligible in B s — ¥ K^~ir + , K^~ p + decays, hence the LL and RR Mis can 
not explain the measurements of B s — > K ~7r + decay from CDF within 2<r ranges. For the LR 
and RL Mis, we have fairly constrained the LR and RL MI parameters from B s — > K~tt + decay 
and B° — B° mixing. Furthermore, using the survived parameter spaces, we have explored the 
LR and RL MI effects on the observables of three B s — > K*~n + ,K ( -*' I ~ p + decays, which have 
not been measured yet. 

The LR and RL insertions can generate sizable effects in B s — > K^~7T + , K^~ p + decays 
since their contributions are enhanced by m-gjm^. The allowed regions of both the SUSY 
weak phases 4>lr,rl and the moduli \(Sf jRRL )i 3 \ have been strongly constrained from B s — > 
K~tt + decay and B° — B° mixing within 2a ranges. We have found as long as x < 2, there 
still is allowed spaces for the LR and RL MI parameters. For x = 0.25 case, B s — > K~n + 
decay and AM^ provide the most stringent constraint, however, for x = 1 case, B s — > K~ir + 
decay and sin2/3 provide the most stringent limit. The theoretical predictions including the 
constrained LR and RL MI contributions are compatible with the measurements within 2cr 
ranges from CDF collaboration in B s — ¥ K~ir + decay. We have found the constrained LR and 
RL insertions still have obvious effects on B(B S -> A'*~7r + , K ( -*^p + ), Acp r (B s -)■ K*~ p + ) and 
fi{B s — > K*~ p + ), moreover, the LR insertion has obvious effects on Aq P (B s — > K~ p + ). Then 
we have presented the sensitivities of the physical observable quantities to the constrained 
LR and RL parameter spaces in Figs. [SB We have found A^ P (B S — > K~p + ,K*~p + ) and 
fh{B s — > K*~ p + ) are very sensitive to the weak phases of the LR and RL insertion parameters. 
In addition, we also have shown the correlations between two observables of B s — > K~ir + decay 
and the observables of B s — ¥ K*~tt + , K^~ p + decays in Figs. EHHJ And we have found all 
observables of B s — > K*~ir + , K^~ p + decays are not very sensitive to B(B S — > K~n + ) for the 
LR and RL Mis, B(B S -> K*~ti + , K^~ p + ) and A% r P (B s ->• K*~p+) have some sensitivities 
to Af P (B s -» K~7i + ) in the case of the LR and RL Mis, and A d £ r P (B s ->■ K*~n + , K~ p + ) are 
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also sensitive to Aqp(B s — > K~ix + ) for the LR MI case. The future measurements or precise 
measurements of the branching ratios, the direct CP asymmetries and the polarization fractions 
in B s — > K^~7r + , K^*^ p + decays could be used to shrink/reveal/rule out the relevant LR and 
RL MI parameter spaces. The results in this paper could be useful for probing SUSY effects 
and searching direct SUSY signals at Tevatron and LHC in the near future. 
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